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Abstract 

We study an ensemble of closed random paths, embedded in M^, with a curvature 
dependent action. Previous analytical results indicate that there is no crumpling 
transition for any finite value of the curvature coupling. Nevertheless, in a high 
statistics numerical simulation, we observe two different regimes for the specific heat 
separated by a rather smooth structure. The analysis of this fact warns us about 
the difficulties in the interpretation of numerical results obtained in cases where 
theoretical results are absent and a high statistics simulation is unreachable. This 
may be the case of random surfaces. 



1 Introduction 

It is well known that the addition of an extrinsic curvature term to the Nambu-Goto action 
{i.e. the area action) for the case of random surfaces controls the formation of spikes, i.e. 
deformations that are exact zero-modes of the area action and causes the degeneration of 
the random surface into branched polymers ||, ||. 

The presence of a second order phase transition for a finite value of the curvature- 
coupling separating a crumpled and a flat phase for the case of fixed connectivity random 
surfaces in three dimensions has been firmly established by recent large scale numerical 
simulations ^ ||. 

Nevertheless, the situation is far from being so clear for the case of dynamically triangu- 
lated random surfaces with extrinsic curvature (DTRS). Results from numerical simulations 
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of this model are unclear and even the real existence of a phase transition separating the 
two regimes has been questioned |^. 

It has been pointed out also that these approaches, initially developed to give a lattice 
discretization for the string theories P], may have a realization in condensed matter physics 
as models of membranes or interfaces. In this sense, fixed connectivity random surfaces are 
called crystalline surfaces and dynamically triangulated surfaces are called fluid surfaces. 
The comparison between the behavior of DTRS and the true fluid membranes states an 
interesting question. In fact, the Monte Carlo analysis of the scaling properties of self- 
avoiding fluid membranes with an extrinsic bending rigidity indicates that fluid membranes 
are always crumpled at sufficiently long length scales ||^, [10[- The observed peak in the 



specific heat is interpreted as the effect originated when the persistence length describing 
the exponential decay of the normal-normal two point function in the crumpled regime 
simply reaches the finite size of the system at the apparent transition. Can this behavior 
be translated to the case of DTRS? 

In this letter we propose a simple scenario: the numerical simulation of the behavior 
of a set of random paths with an extrinsic curvature term in the action. This system is is 
particularly adequate to test the numerical techniques used in random surfaces for several 
reasons: 

1. Its numerical simplicity (not triviality!) that allows to perform simulations with 
several millions of sweeps (simply impossible to do in random surfaces). 

2. It is the most natural place to start with random surfaces simulations, as it has been 
pointed out by several authors [|ll|, [l^, [T^ . 

3. It exist a solid framework given by the analytical results obtained by Ambj0rn et 
al. [Pll. 



2 Random Paths with Curvature 

The analysis of a random path consisting on straight segments embedded in was per- 
formed as a one-dimensional analogue to check the numerical techniques applied in the 
numerical simulations of triangulated random surfaces Pisarski |[T^ proposed the ad- 



dition of an extra term -proportional to the curvature- to the theory of random paths, 
claiming that this theory might be relevant to polymer physics. He noticed the observation 
of asymptotic freedom and pointed out the similarities between this theory and a nonlin- 



ear cr-model with long-range interactions. Furthermore, Alonso and Espriu [13] included 
random paths as a simple version of the random surface theory in a mean field analy- 
sis. They concluded that the Hausdorff dimension of the paths will change from dn = 2 
(Brownian paths) to dn = 1 at infinite curvature coupling. In 1987 Ambj0rn, Durhuus and 
Johnsson [|T^ performed a complete analytical analysis of this model and they concluded 
that for finite values of the coupling constant of the curvature term, the theory belongs 
to the same universality class of simple random walks. A non-trivial scaling limit is only 
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possible if the curvature coupling tends to infinity. They computed exactly the two-point 
function in this limit. They observed that all the paths have anomalous dimension r] = 1, 
which is consistent with the mean field analysis. 

In this letter we perform an actual simulation of a micro-canonical set of random paths 
with a curvature dependent action. Such an analysis, dealing with finite paths made of a 
given number of straight segments, allows for a graphical visualization of the role of the 
curvature term of the action. Geometrical magnitudes such as the curvature or the gyration 
radius can be directly measured and the scaling behavior explicitly deduced. A close 
correspondence with the conventions and techniques used in the simulations of random 
surfaces has been imposed: closed paths of a given number of points, fixed topology, a 
Metropolis algorithm, etc. 



3 Lattice action for closed paths 



Pisarski [|T^] considered a one dimensional analogue of the string action defined in 
allowing a path integration over paths as 



S = l3 J ds + X J k ds, 



(1) 



where s is the parameter for which = 1, and k is the curvature given by 

(both P and x being positive) .[] 

Given an arbitrary parameterization 



s = x 



X2 



dt + (5 iX')^dt, 



(2) 



Pisarski showed, using perturbation theory, that at x^^ = the theory is asymptotically 
free. Therefore, the paths will become smoother as we go to shorter distances, but longer 
distances are outside the validity of perturbation theory. 

Ambj0rn et al. studied the non-perturbative regime of this theory but considered 



open endpoints. Using a lattice-like regularization and a generalized central limit theorem, 
they concluded that there are actually no other phase transitions than the x~^ = point 
(which it is expected from asymptotic freedom). This implies that the correlation length 
must increase for small x^^- A similar result was obtained in the approximation of mean- 
field theory 

The actual action we have simulated is a lattice transcription of the action (P 
by simply replacing the derivatives of the path Xf by finite differences X-^^^ 



obtained 
Xf (and 



doing analogously with second derivatives): 

N 



TV 



s 



=1 



=1 



9,. 



sm 



(3) 



^ As pointed out in the curvature term must be proportional to k, which is scale invariant, and 
requires a dimensionless coupling. Terms proportional to higher powers are formally irrelevant in the 
continuum limit. 
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where 9i is the angle between two adjacent segments. 

The partition function of this theory is, then, defined as 

„ N-l 

Z = j W d'X, C'. (4) 

i=l 

It is easy to check that, hke in the surface case, this action is invariant under reparam- 
eterizations which may be thought as a kind of gauge symmetry. This fact allows to fix 
the coupling j3 = 1 obtaining, actually, a single parameter action. 

We have restricted ourselves to the case of closed paths, i.e. we have require that the 
neighbor of the i = N point is i = 1. We have also fixed this point to avoid translational 
invariance. All these prescriptions are in perfect agreement with the simulation of random 
surfaces. 

From a geometrical point of view, the paths we have simulated consist in closed loops 
of straight segments embedded in M^, being this topology preserved during the simulation. 



4 Numerical results 

We have performed numerical simulations for paths embedded in {D = 3) with different 
number of points: = 25, 100, 200 and for different values of the curvature coupling x 
from to 10, keeping always (3 = 1. The number of sweeps has been of the order of 
2-4 millions for each coupling. The numerical computation of the partition function has 



been done using a simple Metropolis algorithm [|TT] . We have not applied histogram or re- 
summation techniques due to the large fluctuations and large correlations experienced by 
this system. The use of other methods such as Hybrid Monte Carlo algorithms seems not 



to improve the results [|T5] . Main results concerning the geometrical and thermodynamical 



behavior of the paths are as follows: 

1. Length. The measure of the mean length of the paths constitutes a good test of 
the numerical procedure. It is easy to deduce that independently of the curvature x 
value. 

Results to compare with the analytical prediction are summarized in table |I|. 

2. Curvature. The mean total curvature of the path is given by 

We have summarized in fig. ^ the results concerning the curvature measurements. 
Notice the decrease of the mean curvature as the coupling x grows. The solid curve 
represents the limit % ^ oo and N ^ oo when the path becomes a flat circle and 
the total curvature is just 27r. 
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Figure 1: Curvature for several paths 
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Haussdorf dimension (|3=1.0) 
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Figure 2: Hausdorff dimension for several paths 
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3. Hausdorff Dimension. The square gyration radius is a geometrical magnitude 
defined as 

For zero curvature coupling the mean value of the square gyration radius is just [jl6| 



12/?2 ' 



c 



in perfect agreement with our results. In the infinite curvature limit |15 



From the square gyration radius we can extract the Hausdorff dimension (see fig. |^). 

From the analytical results [|T^ we know that in the continuous limit random paths 
must be crumpled for all finite curvature couplings (including the case x = 0). This 
implies dn = 2 (i.e. the square gyration radius scales as A^). Rigid paths are expected 
only at the infinite curvature coupling, and in this case, we expect rf^ = 1 (i.e. the 
square gyration radius scales as A^^). The analysis of fig. ^ shows an unexpected 
behavior. For x = the values of the Hausdorff dimension, for the different sizes 
of the paths, are all higher than dn = 2. These values decrease when increases 
and they are in perfect agreement with eq. so in the infinite volume (A^ oo) 
the value of the Hausdorff dimension is dn = 2, as expected for Brownian paths. 
However, the behavior of paths at non-zero curvature coupling is rather difficult to 
connect with the above scenario. We observe that for each size A^, the value of dn 
decreases as a function of x up to an almost constant value for large (but still finite) 
couplings. For all the sizes studied, these limiting values are lower than dn = 2 and, 
for large A^, they approaches a value close to dn = 1.5. The effects of dealing with 
paths of finite size seems to be, then, stronger than naive expectations. 

We have drawn in fig. ^ simple snapshots of the paths for different curvature cou- 
plings. They show an apparent change of behavior from crumpled or Brownian paths 
to rigid paths. For large curvature terms, a closed path build of a finite number of 
straight segments seems to be frozen. 



Table 1: Some lengths compared with the predicted value 



X 


A^ = 25 (72) 


A^ = 100 (297) 


A^ = 200 (597) 


0.0 


71.9 ±0.1 


296.7 ±0.2 


597.3 ±0.4 


1.4 


71.7 ±0.1 


296.5 ±0.2 


595.9 ±0.5 


3.4 


72.3 ±0.2 


297.5 ±0.6 


600.1 ± 1.0 


5.4 


71.6 ±0.2 


297.5 ±0.8 


594.8 ± 1.6 
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Figure 3: Simple snapshots for several paths. 



4. Specific heat. In addition to this geometrical picture, we have measured also the 
heat capacity with respect to curvature, defined as 



InZ, 



(10) 



where Z is the canonical partition function of the system. 

Its behavior, shown in fig. ^, is really suggestive. Remark, for instance, the saturation 
for long paths. The shape of this graph is very similar to that corresponding to the 
self-avoiding fluid membranes 0, for which the authors conclude the absence of 
a phase transition. Moreover it is also similar to the case of crystalline random 
surfaces, where there is a true phase transition [^, This reflects the fact that a 
crude observation of a structure in the specific heat graph is not relevant enough to 
characterize a phase transition. 

We should remark that a great difficulty in this analysis is the large correlation found 
in the simulations, specially for large curvature couplings. This is not only due to the poor 
efficiency of the method, but is mainly due to the finite size effects of the system. 
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Figure 4: Specific heat for several patlis. 
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5 Conclusions 

We liave performed a numerical simulation of an ensemble of closed random paths em- 
bedded in weighted with an action proportional to both, the length of the path and 
the curvature. Several geometrical and thermodynamical magnitudes have been measured: 
Length, curvature, gyration radius, Hausdorff dimension and specific heat. 

This is a privileged framework, because from analytical results we know that there is 
no true phase transition and that for a finite non-zero curvature coupling these paths will 
belong to the same universality class than Brownian random walks. In addition, one has 
the possibility of perform a high statistics simulation. Nevertheless, in an actual numerical 
simulation, the finite size effects and the large correlations in the system hide the expected 
behavior. 

In addition to the intrinsic interest of such a simulation -that completes previous works 



on this field |12, 14]- the experience acquired warns us that one must be very careful when 



dealing with other situations without firmly established analytical results. 
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